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ABSTRACT 

We present high resolution (0'.'3 = 16 AU), high signal-to- noise ratio Submillimeter 
Array observations of the 870 nm (345 GHz) continuum and CO J=3— 2 line emission 
from the protoplanetary disk around TW Hya. Using continuum and line radiative 
transfer calculations, those data and the multiwavelength spectral energy distribution 
are analyzed together in the context of simple two-dimensional parametric disk 
structure models. Under the assumptions of a radially invariant dust population 
and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously 
reproduce the CO and dust observations with model structures that employ either a 
single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, 
we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a 
relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared 
scattered light) that extends to a much larger radius of at least 215 AU. We discuss 
some possible explanations for the observed radial distribution of millimeter-sized 
dust grains and the apparent CO-dust size discrepancy, and suggest that they may 
be hallmarks of substructure in the dust disk or natural signatures of the growth and 
radial drift of solids that might be expected for disks around older pre-main sequence 
stars like TW Hya. 
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Introduction 



The physical conditions of the gas and dust in young circumstellar disks shape the formation 
and early evolution of planetary systems. The spatial distribution of disk material - the density 
structure ~ plays many fundamental roles, and is especially relevant for dictating when and where 
planets can form and migrate in the disk. Naturally, empirical constraints on protoplanetary disk 
structures are of significant value in efforts to develop realistic models of the complex processes 
involved in planet formation. Some recent studies have made progress in extracting the radial (sur- 
face) density profiles in the se disks from high arigular resolution measurements of their millimeter - 



wave continuum emission (Andrews et al. 200S 



2ninbl : llsella et al.l 120091 : lOuilloteau et al.l I2OII 



However, such observations are only sensitive to the density structure of the (presumably) trace 
population of disk solids, and not the gas that wields far more influence over the evolution of disk 
structure and the formation of planets. The major challenge is that the gas in these disks is pri- 
marily cold H2, which is not directly observable. And wh ile there have been many investigations 
of less abundant molecular species (see lDutrev et al.l 120071 ). limited sensitivity to low optical depth 
lines and an incomplete understanding of the complicated chemistry in th ese disks have so far mad e 



it difficult to convert observations of the gas into density constraints (see IWilliams &: Cieza 



2011 



Given those obstacles, there have not been many attempts to directly compare the gas and dust 
structures of protoplanetary disks with spatially resolved data. A few studies used simple models to 



(Pietu et al. 


2005; 


Isella et al. 


2007) 



Hughes et al.l (|2008l ) suggested that this discrepancy was likely an optical depth illusion, an artifact 
of the assumed sharp outer edge in the density model. They showed that a surface density profile 
with a smoo t h tap er at large radii could better reproduce the dust and gas observations. However, 



Panic et al 



(120091^ ar gued that similar modifications to the model structure of the IM Lup disk 



(see iPinte et al.l 120081 ) were insufficier it to account fo r the much different sizes they measure for 
its gas and dust emission. Recently, iQi et al.l (|201ll ) successfully matched their observations of 
multiple CO transitions and continuum emission from the HD 163296 disk with a single density 
model, although refinements to that model based on the detailed morphology of the continuum 
emission were not a priority in that study. In general, there is still substantial uncertainty that the 
gas and dust trace the same structures in protoplanetary disks. From a practical standpoint, that 
uncertainty is disconcerting. Our knowledge of the mass contents of these disks is entirely based 
on their (easy to measure) dust emission, but we are not sure if those dust structures also describe 
the gas reservoirs that effectively control the key disk evolution and planet formation processes. 

Although direct measurements of H2 densities are not possible, in principle an indirect com- 
parison of the dust and gas structures in a protoplanetary disk can be made with sufficiently 
sensitive and high resolution observations of the thermal continuum and emission lines from trace 
gas species. To meet those requirements, a relatively massive and nearby disk would be an ideal 
candidat e target. TW Hya is an i s olated, ^0.8 Mf;^ T Tauri star located on ly 54 it 6 pc from 
the Sun (jRucinski fc Krautterl Il983l : IWichmann et al.l Il998l : Ivan LeeuwenI 120071 ) . Despite its ad- 
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vanced age (^lOMvr: iKastner et al.lll99?l : IWebb et al.lll999l ). it hosts a massive, gaseous accre- 



tion disk with a rich molecular spectrum and strong cq n tinuu n i emi ssion out to centimeter wave- 



lengths (jWilner et alJ l200d . l2003l . l2005l : lOi et alJ 12004 . l2006l . 120081 ) . Those observations and a 



viewed nearly face-on ( 


Krist et al. 


2000: 


2004 




Roberge et al. 


2005"). The inner ed 



4AU, likely by an unseen (and maybe planetary) companion (jCalvet et al 



2002 



Hughes et al 



20071 ). Given its proximity, favorable orientation, and rich gas and dust content, the TW Hya disk 
is uniquely well-suited for a comparative investigation of the radial behavior of its gas and dust 
structures. 

In this article, we present new sub-arcsecond resolution observations of the 870 fim continuum 
and CO J=3— 2 emission from the TW Hya disk. Using state-of-the-art radiative transfer modeling 
tools, we use those measurements to compare the radial distributions of its CO gas and dust. In §2 
we describe the observations and data calibration. The basic characteristics of the data are reviewed 
in §3. A detailed description of the modeling and results are provided in §4. The implications of 
that analysis are discussed in §5, and our principle conclusions are summarized in §6. 



Observations and Data Reduction 



TW Hya was observed at 345 GHz (870 /^m) with the Submillimeter Array (SMA; iHo et al 



20041 ) on 8 occasions in 2006-2010 in the compact (C), extended (E), and very extended (V) array 
configurations, providing baseline lengths of 16-70 m, 28-226 m, and 68-509 m, respectively. To 
accomodate various studies of the CO J=3— 2 emission line (at 345.796GHz), the SMA double 
sideband receivers and correlator were configured with several different local oscillator (LO) settings 
and spectral resolutions throughout this campaign. The "default" setup used an LO frequency of 
340.755 GHz (880 fim) with the CO line in the center of the upper sideband, with a resolution of 
0.70 km s~^ in one spectral chunk. The other 23 partially-overlapping 104 MHz spectral chunks 
in each sideband were sampled coarsely with 32 channels each. Some of the 2008 observations 
utilized a "high" resolution correlator mode, with a large portion of the bandwidth devoted to 



probing the CO line with 0.044 km s ^ channels (see Hughes et al.|[201ll ). While those data sample 



the continuum with the default spectral resolution, they have a reduced continuum bandwidth 
(2.6 GHz; ~70% of the available bandwidth at the time) and a higher LO frequency at 349.935 GHz 
(857 /um). The 2010 observations were conducted in a "medium" resolution mode that employed the 
new expanded bandwidth capabilities at the SMA. In that case, two different 2 GHz IF bands were 
centered ±5 GHz (the normal setup) and ±7 GHz from the LO frequency, 339.853 GHz (883 //m). 
The CO line was observed at a moderate resolution of 0.18 km s^^ in the upper sideband of the 
first IF band, and the continuum was coarsely sampled across the remaining expanded bandwidth. 
Fin a,lly, a "hybrid" mode was utilized at the end of 2006 to cover the J=4-3 transition of H^3qo+ 



(see lOi et al.ll2008l ). A summary of the observational setups is provided in Table [TJ 
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Table 1. Summary of SMA Observations 



UT Date 


array 


spectral 


RMS noise 


beam size 


beam PA 




config. 


setup 


[mJy beam" 


'] ["] 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2006 Dec 28 


C 


hybrid 


3.7 


4.2 X 1.7 


4 


2008 Jan 23 


E 


high 


3.2 


0.9 X 0.7 


21 


2008 Feb 20 


E 


high 


3.1 


1.0 X 0.7 


-20 


2008 Feb 21 


E 


default 


4.1 


1.2 X 0.5 


14 


2008 Mar 2 


C 


high 


3.6 


3.9 X 1.9 


-17 


2008 Mar 9 


C 


default 


3.2 


3.7 X 1.8 


-13 


2008 Apr 3 


V 


default 


2.2 


0.6 X 0.3 


19 


2010 Feb 9 


V 


medium 


1.8 


0.5 X 0.3 


-3 


combined 


all 




2.0 


0.8 X 0.6 


-3 


Note. — Col. (1): UT date of observations. 


Col. (2): SMA 


array con- 


figuration, C : 


= compact 


, E = extended, and V 


= very extended. Col. (3): 


Adopted correlator mode (see §2' 


). Col. (4): Continuum RMS 


noise in a 


naturally-weig 


;htcd map. 


Col. (5): 


Naturally-wei; 


ghted synthesized beam di- 



mensions. Col. (6): Naturally- weighted synthesized beam orientation (ma- 
jor axis position angle, measured east of north). 
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TW Hya was observed in an alternating sequence with the nearby quasar J1037-295 (a pro- 
jected separation of 7.3°), with a total cycle time of 15 minutes in the C and E configurations and 
8-10 minutes in the V configuration. The quasars 3C 279 or J1146-289 were observed every other 
cycle. For the tracks on 2008 February 21, March 9, and April 3 , the target portion o f the cycle 



was shared with the nearby sources HD 98800 and Hen 3-600 (see [Andrews et al.ll2010al ). Planets, 
satellites, and bright quasars were observed as bandpass and absolute flux calibrators when TW 
Hya was at low elevations (<18°), depending on their availability and the array configuration. Most 
of these data were obtained in the best atmospheric conditions available at Mauna Kea, with zenith 
opacities of 0.03-0.05 at 225 GHz (0.6-1.0 mm of precipitable water vapor) and well-behaved phase 
variations on timescales longer than the calibration cycle. The conditions on 2008 March 2 and 9 
were worse, but typical for Mauna Kea (with 1.4-1.8 mm of precipitable water vapor). 

The data from each individual observation were reduced independently with the IDL-based MIR 
software package. The bandpass response was calibrated with observations of bright planets and 
quasars, and broadband continuum channels in each sideband (and IF band, where applicable) were 
generated by averaging and then combining the central 82 MHz in each line-free spectral chunk. 
The visibility amplitude scale was derived from observations of Uranus, Titan, Callisto, or Vesta, 
with a typical systematic uncertainty of 10-15%. The antenna-based complex gain response of 
the system was determined with reference to J1037-295. The observations of 3C 279 or J1146-289 
were used to assess the quality of phase transfer in the gain calibration process. Based on those 
data, we estimate that the "seeing" generated by atmospheric phase noise and any small baseline 
errors is small, 0.10-0'.'l5. This result is a testament to the high quality of the observing conditions, 
especially given the low target elevation and wide projected separations between the calibrators (for 
reference, 3C 279 is located 39.7° from TW Hya and 38.5° from J1037-295, while the corresponding 
separations for J1146-289 are 11.0° and 15.1°, respectively). 

Before the observations could be combined, we had to account for the proper n iotion of TW 
Hya over the ~3 year observing baseline {fiaCosS = — 0'.'066 yr~^, fis = — 0^.^014 vr~^: Ivan Leeuwen 



20071 ). Fortunately, each individual dataset exhibited bright, symmetric, centrally-peaked contin- 
uum emission that we associate with dust in the TW Hya disk. Using elliptical Gaussian fits to 
those continuum visibilities, we determined centroid positions for each individual set of observations 
and associated them with the stellar position at that epoch. The measured emission centroids are 
consistent with the expected stellar positions, well within the ~0'.'l absolute astrometric accuracy 
of the SMA (set primarily by small baseline uncertainties). Based on those centroid measurements, 
each individual dataset was aligned to a common coordinate system. After confirming that the 
aligned visibility sets showed excellent agreement on all overlapping baselines, all datasets were 
combined to produce composite 345 GHz continuum and CO J=3— 2 visibilities. 

The composite visibilities were Fourier inverted, deconvolved with the CLEAN algorithm, and 
restored with a synthesized beam using the MIRIAD software package. The natural weighting of the 
continuum data produces an image with a 0'.'80 x 0'.'58 beam and an RMS noise level of 2.0 mJy 
beam~^ (which is dynamic-range limited). Composite CO J=3— 2 channel maps were synthesized 
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with a velocity resolution of 0.2 km s ^ and a circular beam with a FWHM of I'.'O (the naturally- 
weighted resolution is C/OG x 0'.'73). The RMS noise level is 0.13 Jy beam~^ in each channel. 

Ancillary information in the literature was used to construct the TW Hya spectral energy 
di stribution (SEP) that will be used with the SMA data. We adopted the optical m onitoring results 
of Mekkaden ( 1998 ) and near-infrared photometry from iWeintraub et al. (2000) and the 2MASS 
po int source catalog (ICutri et al.l 120031). I n the thermal infrared, Spitzer flux densities measured 
bvlHartmann et al.l (120051) andlLow et al.l (|2005l ) were supplemented with IRAS and Herschel data 



( Weaver &: Jones 



1992 



Thi et al 



2010l ). We also include a Spitzer IRS spectrum, kindly provided 



by E. Furlan (see lUchida et al.ll2004l ). At submillimeter wavelengths, we relied on the calibrator flux 
densities from the SHARC-lf] and SCUBaI instruments at the Caltech Submillimeter Observatory 
and James Clerk Maxwell Telescope, respectively {F^, = 6.13 ± 0.68, 3.9 ± 0.7, and 1.37 ± 0.01 Jy 
at 35 0, 443, and 869 A ^ m, re s pectively) . Add itiona l millimeter-wave rn easurements were collected 
from lWeintraub et al.l (|l989l ). lQi etaP (|2004l ). and lwilner et all (|2003l ). 



3. Results 

The astrometrically aligned and combined SMA data are shown together in Figure [H along 
with the broadband SED. The synthesized continuum image in Figure [T]2 has an effective frequency 
of 344.4 GHz (870 /um), with an integrated flux density of 1.34 it 0.13 Jy and a peak flux density of 
0.337ib 0.034 Jy beam~^ (a peak signal-to-noise ratio of ~170), including the systematic calibration 
uncertainties. The continuum emission is regular and symmetric on the angular scales probed here, 
with a synthesized beam size of 43 x 31 AU projected on the sky. Essentially no emission is detected 
outside a ~1" radius. Given that the integrated flux de nsity determined from the SMA data is in 



good agreement with single-dish measurements (see §2: IWeintraub et al.lll989l : iDi Francesco et al 
20081 ). it is clear that this is no optical depth effect: all of the millimeter- wave dust emission from 
the TW Hya disk is concentrated inside a projected radius of ~60 AU. 

While the continuum image appears rather plain, some interesting emission features are appar- 
ent in a direct examination of the visibilities. Figure [1^ displays the azimut hally averaged 870 um 
continuum real visibilities as a function of the deprojected baseline length (see I Andrews et al.ll2009l ). 
assuming the inclination (i = 6°) and major axis position angle (PA = 335°) derived from high 



spectral resolution CO emission line data by iHughes et al.l (|201ll ). The imaginary visibilities are 



effectively zero within the noise on all sampled baselines, confirming the accuracy of the astromet- 
ric alignment and reinforcing that there are no obvious departures from axisymmetry in the TW 
Hya disk. The visibility profile in Figure [!][) shows a smooth decrease out to ~150kA, followed by 
low-amplitude oscillations on longer baselines and an apparent null near 500 kA. These features are 



' [http: //www, submm. caltech. edu/$\sim$sharc/analy sis/calibration. htm| 

^http: //www, jach.hawaii . edu/JCMT/continuvun/calibratioii/sens/poteiitialcalibrators .html| 
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Fig. 1. — (a) Naturally weighted composite image of the 870 fim continuum emission from the TW 
Hya disk. Contours start at 10 mJy beam""*^ (5 a) and increase in 30 mJy beam~^ (15 a) increments. 
The synthesized beam dimensions are shown in the lower left. (6) Azimuthally averaged 870 /um 
continuum visibility profile as a function of the deprojected baseline length (real part only; the 
imaginary terms are effectively zero on all baselines). The uncertainties are typically smaller than 
the symbol sizes. Note the low-amplitude oscillations beyond ~150kA. (c) Complete SED for 
the TW Hya star-l-disk system (references arc in the text; sec §2). The Spitzer IRS spectrum is 
shown as a thick gray curve. Our adopted stellar photosphere model is overlaid as a thin gray 
curve (see §4.1). (d) CO J=3— 2 channel maps from the TW Hya disk, on the same angular scale 
as the continuum map in panel a. Contours are drawn at 0.4 Jy beam~^ intervals (~3(j) in each 
0.2 km s~^-wide channel. The 1" synthesized beam is shown in the lower left. The central channel 
represents the TW Hya systemic velocity, at Vlsr = +2.86 km s"^. 
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distinct in independent datasets (particularly the "dip" near 180 kA, where 4 different observations 
span that range of baselines), are present regardless of the bin sizes used for profile averaging, 
and are not noted in the visibility profiles for the test calibrators (J1146-289 or 3C 279). Despite 
the challenges of calibrating SMA data for low-elevation targets, the persistence of these visibility 
modulations make us confident that they are real features. Moreover, we will demonstrate in §4 
that they can be reproduced with models that incorporate a sharp outer edge in their emission 
profiles. The null at 500 kA is consiste nt with the 4AU-radius inner disk ca vity inferred from t he 
SED (see Figure [Tb: ICalvet et al.ll2002l ) and a VLA 7 mm continuum image (jHughes et al.l 120071 ). 



The panels in Figure[T]i show the composite CO J=3— 2 emission line channel maps for the TW 
Hya disk, resampled to a velocity resolution of 0.2 km with a circular I'.'O (54 AU) synthesized 
beam. The emission is firmly detected (>3o") out to ±1.2 km s~^ from the systemic velocity 
(^LSR = 2.86km s~^), with an integrated intensity of 34.8 it 3.5 Jy km s~^ and a peak flux of 
4.2 zb 0.4 Jy beam~^ (43 ± 4K), including the calibration unce rtainties. Those values are in good 



agreement with previous single-dish and SMA measurements (jvan Zadelhoff et al 



2004 



Hughes et al 



2011 



2001 



Qi et al 



The channel maps in Figure [T]i show a clear rotation pattern, from 
northwest (blueshifted) to southeast (redshifted), with a narrow line- width due to a face-on viewing 
geometry. Near the systemic velocity, the CO emission subtends ~4" (215 AU) in radius. 



4. Modeling Analysis 



These SMA observations offer some new insights into the TW Hya disk structure. Nat- 
urally, as one of the nearest pre-main sequence stars, TW Hya and its associated disk have 
been the subject of intense observational scrutin y. The global stru cture of the TW Hya disk 



Roberge et a 



Hughes et al 



2008 



2011 



light observations (|Krist et al.l l2000l : iTrilling et al.l 1200 ll : IWeinberger et al 



200a), millimeter /radio- wave continuum images (IWilner et al 



2OO7I 1. and resolved molecular line maps (jOi et al 



2002 



: Aoai et al. 


200J; 


2OO0I. 


2OO3I, 


2OO5I; 



2OO4I . bood . I2OO8I : iHuehes et al" 



). However, none of those previous studies had the combination of angular resolution 
and sensitivity for the optically thin dust and high-quality gas tracers that are available from the 
SMA datasets presented here. 

In the following, a technique is described for extracting the structure of the TW Hya disk from 
these data using radiative transfer models. We focus specifically on enabling a comparison between 
the radial distributions of the dust and CO tracers. The approach we have adopted to make that 
comparison consists of three key steps. First, we construct a model of the radial density structure 
of the dust that is able to reproduce our resolved observations of 870 fim continuum emission and 
the broadband SED. Next, we make the assumption of a radially constant (vertically-integrated) 
dust-to-gas mass ratio and use the model structure we derive from the dust to predict the emission 
morphology of the CO J=3— 2 line. Then, we show that this assumption implies an inconsistency 
with the observations, highlighting a clear difference in the radial distributions of millimeter-sized 



- 9 - 



dust grains and CO gas in the TW Hya disk. Some of the potential impUcations of this inconsistency 
are discussed further in S5. 



4.1. Dust Structure 



The dust disk structure is determined following the technique outlined by [Andrews et al 



(|201ll ). with some modifications for generality. We assume the dust is spatially distributed with a 
parametric two-dimensional density structure in cylindrical-polar coordinates {r, z}, 



Pd{r, z) 



2'KZd 



exp 



(1) 



where and z^ are surface densities and characteric heights, which both vary radially (see below). 
As will be explained further in §4.3, we investigated two different models for the radial surface 
densit y profile. First, we employed the similarity solution for simple viscous accretion disk struc- 
tures (|Lvnden-Bell &: Pringldll974l) that we have used su cc essfully to characterize bo th normal and 
transition disks in the past (jAndrews et allbood . boiOal lbl. [ 201ll : iHughes et al.ll20ld ). In that case. 



exp 



2-7 



(2) 



where Sc is a normalization, Vc is a characteristic scaling radius, and 7 is a gradient parameter. 
As an alternative, we considered a less physically motivated (but perhaps more commonly used ) 
model that incorporates a power-law density profile with a sharp cut-off (see lAndrews et al.ll2008l ). 



Sd(r) 



(if r < ro; else = 0), 



(3) 



where Sq is a normalization, rg is the outer edge of the disk, and p is a gradient parameter. In 



(Calvet et al. 


2002; 


Hughes et al. 


2007) 



we set the surface densities to a constant value Sin between the sublimation radius (rsub) and a 
"gap" radius (^gap). No dust is present between that gap radius and the cavity edge, rcav In the 
vertical dimension, the dust is distributed like a Gaussian w ith a variance z^. Th e characteristic 

we employ a 



(201: 



height varies with radius like z^ = ZQ{r /ro)^^^ . Following [Andrews et al. 
cavity "wall" to reproduce the infrared spectrum of TW Hya (no such feature was required at the 
sublimation radius). The local value of z^ is scaled up to z^aii at rcav, and then exponentially joined 
to the global z^ distribution over a small radial width, Ar„aii- 

This structure model has 11 parameters: three describe the base surface density profile, {Sc, 
rc, 7} or {So, ro, p}, five determine the cavity and inner disk properties, {Sjn, rsub, ^gap, ^'cav, 
Arwaii}, and three others characterize the vertical distribution of dust, {zq, z^^w-, '4'}- To simplify 
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the modeling, we fixed some of the parameters that are of less direct interest here. The sublimation 
radius was set to rgub = 0.05 AU, the location where dust temperatures reach 1400 K (see also 
Eisner et al.ll2006l ). The gap radius was set to ^gap — 0.3 AU and the (constant) inn er disk density 



to Ein = 5x 10 ^ g cm ^. The cavity edge was fixed at rcav = 4 AU (see lHughes et al.ll2007l ). the wall 



height was set to z^aii = 0.25 AU, and the wall width to Ar- 



wall 



1 AU. Since the details of this gap 



Eisner et al. 


2006; 


Ratzka et al. 


2007 


Akeson et al. 


2011 



with modeling the SED, we also fixed the scale height gradient to ip = 0.25. The in t erpla. y and 



degeneracies between these free parameters were discussed in detail by [Andrews et al.l (|201ll ). For 
our purposes here, it is worth emphasizing that the parameters we have fixed have little quantitative 
impact on the derived radial structures (i.e., sizes and density gradients). 



We used the dust composition advocated by iPollack et al.l (jl994l ). consisting of a mixture of 
astronomical silicates, water ice, troilit e, and organics with th e abundances, optical p roperties, and 



sublim ation temperatures discussed by lD'Alessio et al.l (120011 ). Based on the efforts of lUchida et al 



(j2004l ) to faithfully reproduce the details of the Spitzer IRS spectrum, we let 25% of the total 
silicate abundance ins ide the disk cavit y (r < rcav) be composed of crystalline forsterite (using 
optical constants from lJager et al.ll2003l ). Two grain populations were employed, with a power-law 
size (s) distribution, n(s) oc s~^'^, between Smm = 0.005 /um and a given Smax- Outside the cavity 
wall, 95% of the dust (by mass) has Smax = 1 mm and the remaining 5% has Smax = 1 /^ni- The 
dust in the wall itself and the tenuous inner disk was assumed to have Smax = 1 /^m- No effort was 
made to distinguish the vertical distributions of these dust populations. Opacity spectra for each 
population were determined from Mie calculations, assuming segregated spherical grains. For these 
dust assumptions, the 870 /im dust opacity in the outer disk is K^nm = 3.4 cm^ g~^. 

We assumed the central star has a K7 spectral type with Tpff = 4110 K, .R^ = 1.04 R©, and 
-M* = 0.8 Mq (logs' = 4.3), based on an effort to match iLeieune et al.l (119971 ) spectr al synthesis 
models to the broadband SED and the detailed optical/infrared spectral analysis work of lYang et al 



(|2005l ) . The bes t-fit stellar spectrum te mplate is overlaid on the SED in Figure [T]i as a light gray 
curve. Recently, Vacca &: SandelJ ( 2011 ) have argued instead for a M2.5 spectral type in the near- 
infrared, and a correspondingly cooler stellar photosphere (3400 K), larger radius (1.29 Rq), and 
lower mass (0.4 Mq). While that stellar model provides a good match to the broadband infrared 
photometry for TW Hya, it underpredicts the observed optical fluxes by a factor of ^3 ( in the 
BVR bandpasses, and its known variability does not bridge that gap; see lMekkadenlll998l ). We 
prefer the parameters for the warmer photosphere because they produce a template spectrum that 
better matches the SED across the complete set of optical and infrared bandpasses. 

For a given set of parameters, we simulated the stellar irradiation and emission output of a 
model dust structure using the twq-dimeri sional, axisymmetric Monte Carlo radiative transfer code 
RADMC (see lEjullemond &: Dominild l2004al ) . Assuming the fixed viewing geometry determined by 
Hughes et al.l (120 111 ), a raytracing algorithm was then used to compute a synthetic model SED 
and set of 870 ^m continuum visibilities sampled at the same spatial frequencies observed with the 
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SMA. For each surface density model, we found the best simultaneous fit to the observed SED and 
SMA visibilities over a coarse grid of the gradient parameter 7 or p, by varying the parameters {Ec 
or So, Vc or rg, zq}. Based on those results, we refined our search and permitted the gradients to 
vary freely to find the best-fit parameter sets for each model type (see §4.3 for results). 



4.2. CO Gas Structure 



Unlike for the dust, the radial density profile of the gas disk cannot be inferred directly from 
models of the optically thick J=3— 2 transition of CO. Therefore, we make a fundamental assump- 
tion that the gas traces the dust in the radial dimension. For any given S^, we define the gas 
surface density profile as T,g = S^/C, where ^ is a (radially) constant dust-to-gas mass ratio. 

However, we have elected to permit some freedom in the vertical distribution of the gas to 
facilitate a more faithfu l reproduction of the CO channel maps. Using a multi-transition CO 
dataset, Qi et al. (j2006l ) noted that models of the TW Hya disk structure had a difficult time 
reproducing an appropriate vertical temperature gradient of the gas. The intensity of the high- 
excitation J=6— 5 line indicated that the gas in the disk atmosphere was significantly hotter than 
the dust, presumably due to substantial X-ray heating from the central star. To be able to reproduce 
the observed CO spectral images, we have characterize d the vertical t empe rature profile of the gas 
in parametric form, based on the modeling analysis of lDartois et al.l (j2003l ). We assume that 



Tg{r, z) =Ta + {Tjn - To) cos ( ^ 



2<5 



(4) 



where Ta = Ti(r/1 AU)~^ is a parametric radial temperature profile in the disk atmosphere, is 
the midplane temperature determined from the RADMC simulations of the dust, 5 describes the shape 
of the vertical profile, and Zq defines the height of the atmosphere layer such that Tg[z > Zg) = Ta- 
in our modeling, we fix 5 = 2 and Zq = 4:Hp, where Hp is the pressure scale height assuming the 
midplane temperature at each radius {Hp = c^/O, the ratio of the midplane sound speed to the 
Keplerian angular velocity). In practice, Eq. (4) is a r easonable paranietric approximation of the 
vertical temperature profile in an irradiated disk (e.g., iD'Alessio et al.lll999l ). We have grounded 



the models by forcing Tg = at the midplane {z = 0), but allowed the gas temperatures to 
increase faster with height than the dust to simulate any additional external heating sources. For 
a given Tg{r, z) specified by {Ti, q}, we then calculate the vertical density structure of the gas by 
numerically integrating the equation of vertical hydrostatic equilibrium. 



din 



Pa 



dz 



GM^:Z \ / iimH\ 91nT( 



(r2 + z2)3/2 J \kTg J 



+ 



dz 



(5) 



using Eg as a boundary condition, where G is the gravitational constant, /x = 2.37 is the mean 
molecular weight of the gas, rriH is the mass of a hydrogen atom, and k is the Boltzmann constant. 
For reference, a vertical slice of a representative model structure is shown in Figured 
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Fig. 2. — Schematic demonstration of our model vertical structure, shown as a vertical slice at 
a fixed radius, (top) The density profile of the gas (red) and dust (blue) as a function of height 
above the midplane. The latter has a parametric Gaussian distribution with a variance z'^, while 
the former is computed assuming it is in hydrostatic pressure balance for its specified temperature 
structure. The relative normalizations of each are represented accurately, such that = C'^g, 
where C is a fixed dust-to-gas mass ratio. The hatched region near the midplane {z < Zm) marks 
where CO is depleted from the gas phase because it is frozen onto dust grain mantles (where 
Tg < Tfrz). The comparable surface CO depletion zone due to photodissociation {z > Zg) is well off 
the right of the plot, (bottom) The corresponding temperature profiles, where is computed from 
the RADMC radiative transfer calculations and Tg is determined parametrically, as described in the 
text. The gas and dust temperatures are equivalent (Tg = = Tm) in the midplane, but Tg rises 
more rapidly than the dust before it saturates to a value at a height Zq. 
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To quantify the de nsity of CO molecule s from any given gas model structure, we adopt 
the layered approach of iQi et alj (|2008l . 120111 ). Based on the detailed chemical calculations of 
Aikawa fc: Nomural (j2006l ). we define two vertical boundaries {zm, Zg} at any given radius such that 
the CO mass fraction (relative to H2) is Xco if Zm{f) ^ z < Zs{r) and 10~^Xco elsewhere. The 
"midplane" boundary, z^, marks the maximum height where CO molecules are expected to be 
frozen out of the gas phase and affixed to dust grain mantles. In practice, Zm is defined as the 
minimum height where Tg > Tfr^, where the freeze-out temperature Tfrz is a parameter considered 
to be constant with radius. The "surface" boundary, Zg, is meant to rep resent the heigh t where CO 
molecules can be photodissociated by X-rays or cosmic rays. Following iQi et al.l (|2008l ). we define 
Zg such that 

1 

fJ-n^H Joo 

where A'pd is a vertically-integrated H2 column density that effectively represents the penetration 
depth of the photodissociating radiation field: A^pd is treated as a radially constant parameter. 

So, for any dust model a corresponding CO model can be characterized with a set of six 
additional parameters: two describe the abundances of dust and CO relative to the total gas mass, 
{C, Xco}, two others characterize the gas temperatures in the disk atmosphere, {Ti, q}, and the 
last two define the spatial distribution of CO in the gas phase, {Tfrz, Apd}- For a given set of these 
parameters, WG gGncratc s. two-dimensional grid of tIqq (r, z) and Tg{r, z) values and define a velocity 
field based on Keplerian rotation around a po int mass M^, , assurn ing a minimal turbulent velocity 
line width of 10m s~^ based on the analysis of lHughes et al. ( 2011). We then feed that information 
into the radiative transfer modeling code LIME (|Brinch Hogerheiiddl2010l ) to solve the non-LTE 
molecular excitation conditions of the model and generate a synthetic high-resolution model cube 
for the CO J=3— 2 transition. That model cube was then re-sampled at the velocity resolution of 
the data, and its Fourier transform was sampled at the same spatial frequencies observed by the 
SMA. In practice, we fixed the dust-to-gas r atio based on the a s sumed dust composit ion, where 
C = 0.014 (a gas-to-dust mass ratio of 71, see lPollack et al.lll994l : iD'Alessio et al.ll200lh . For each 
dust model, we varied {Ti, q} for a coarse grid of CO abundance layer parameters {Xco, Tf^z, 
Apd} to find the best available match to the SMA spectral visibilities. Since we are using only a 
single CO transition in this investigation, the abundance layer parameters do not have a strong, 
independent effect on the synthetic CO visibilities. For simplicity, we adopt the best-fit models 
where = 2 x 10^^, Tfj-^ = 20 K, and Np^^ = 10^^ cm~^ as representative. 



4.3. Modeling Results 

The best-fit parameters for a range of representative surface density gradients of each model 
type are compiled in Table [2j For clear notation, each model is labeled with a letter designation 
corresponding to its 7 or p value, preceded by a lower-case 's' for the similarity solution models 
(based on Eq. 2) or 'p' for the power-law models (based on Eq. 3). The surface densities, char- 



-14- 



Table 2. Estimated Model Parameters and Fit Results 



Model 


S^flO) 


1 5 f 


Tn 
' C? ' U 








Ased 


Acont 


A CO 




[g- rin"^] 


















(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


sA 


0.29 


1.0 


35 


0.62 


104 


0.55 


201 


313,914 


319,974 


sB 


0.51 


0.5 


43 


0.58 


101 


0.54 


204 


304,303 


319,989 


sC 


0.28 


0.0 


45 


0.57 


98 


0.54 


210 


302,132 


320,013 


sD 


0.20 


-0.5 


42 


0.53 


103 


0.53 


216 


304,329 


320,024 


sE 


0.14 


-1.0 


40 


0.50 


102 


0.53 


231 


307,785 


320,040 


pA 


0.79 


1.5 


77 


0.62 


104 


0.54 


205 


309,370 


320,056 


pB 


0.46 


1.0 


67 


0.60 


99 


0.54 


211 


302,747 


320,047 


pC 


0.39 


0.75 


60 


0.58 


99 


0.54 


215 


301,785 


320,044 


pD 


0.31 


0.5 


58 


0.58 


100 


0.54 


227 


302,570 


320,035 


pE 


0.23 


0.0 


51 


0.55 


104 


0.54 


249 


306,980 


320,068 



Note. — Col. (1): Model designation, where 's' = similarity solution and 'p' = 
power- law with sharp edge as defined in Eq. 2 and 3, respectively. Col. (2): Dust sur- 
face density at r = 10 AU. Col. (3): Surface density gradient. Col. (4): Characteristic 
scaling radius ('s' models) or outer edge radius ('p' models). Col. (5): Characteris- 
tic dust height at r = 10 AU. Col. (6): Gas temperature in the disk atmosphere at 
r = 10 AU (see Eq. 4). Col. (7): Radial gradient of the atmosphere gas temperature 
profile. Col. (8): statistic for the SED (including Spitzer IRS spectrum). Col. (9): 
statistic for the 870 /xm visibilities. Col. (10): statistic for the CO J=3— 2 
spectral visibilities. The quoted parameter values and fit results arc valid for a set of 
additional fixed parameters, including: the gradient of the dust height profile = 0.25, 
the sublimation radius rg^\~, = 0.05 AU, the inner ra,(iius of the g3<p ^gap — 

0.3 AU, the 

cavity edge rcav = 4AU, the cavity wall height z^aii = 0.25 AU, the gas temperature 
profile parameters S = 2 and Zq = -iHp, the dust-to-gas ratio ( = 0.014, the CO/H2 
abundance ratio Xco = 2 x 10""^, the CO freezeout temperature Tf^z = 20 K, and 
the CO photodissociation column Np^ = 10^^ cm~^. Furthermore, we assume a disk 
inclination of 6°, major axis position angle of 335°, and stellar mass of 0.8 M©. 
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Fig. 3. — {top) Best-fit gas surface densities (where T,g = T,d/C)- Similarity solution models are 
shown on the left and power-law + sharp edge models are shown on the right, (bottom) Temperature 
profiles at the disk midplane (Tm) are shown as solid curves, and were determined from RADMC 
radiative transfer calculations. The parametric atmosphere temperatures (T^; at a height Zg, see 
§4.2) are overlaid as dashed curves. All models have a gap at r = 4AU, marked with a dotted 
gray line. The CO freezeout temperature, Tfj-z = 20 K, is marked as a horizontal gray line. The 
temperature profiles shown here are similar for all models: the midplane values do not change 
much because the total irradiated dust mass is roughly the same in each case, and the atmosphere 
temperatures are determined from the same CO emission data. However, each model does have 
slight variation in the vertical temperature profile, which is not displayed here for the sake of clarity. 
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Fig. 4. — A comparison of the model structmes in Table [2] with observations of the dust continuum 
emission from the TW Hya disk, including the 870 //m visibility profile {top, with a zoomed-in view 
of the emission on longer baselines in the middle) and SED (bottom). The similarity solution surface 
density models (Eq. 2) are shown on the left, and the power-law models with sharp outer edges 
(Eq. 3) are shown on the right. The model SEDs are essentially indistinguishable, but there are 
clear differences in the 870 fim radial emission profiles. The overall best match to the continuum 
emission is Model pC, which has a density gradient p = 0.75 and a sharp outer edge at ro = 60 AU 
(for the similarity solution models, the best match is Model sC). 
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acteristic dust heights, and gas atmosphere temperatures are Usted for a radius of 10 AU, where 
the two surface density model types have comparable behavior. The four parameters describing 
the dust densities, {£^(10), 7 or p, Tc or ro, 2^(10)}, were determined from joint fits to the SED 
and 870 //m continuum visibilities. With those parameters fixed, the gas temperature parameters 
{Ta(lO), q\ were estimated from the CO J=3— 2 visibilities only. The individual values for the 
SED, continuum visibilities, and CO visibilities are listed in Cols. (8)-(10). There are 77 indepen- 
dent datapoints in the SED (including 45 equally-spaced points across the Spitzer IRS spectrum), 
166,796 continuum visibilities, and 83,740 CO visibilities in each of 15 spectral channels (a total 
of 1,256,100). The values were calculated with weights that incorporate the quadrature sum 
of the formal uncertainties and absolute calibration uncertainties on each datapoint (e.g., a 10% 
systematic uncertainty on the amplitudes) . Figure [3] shows the radial structures for each of the 
disk models in Tabled including their surface densities {top) and temperature profiles (bottom). 

Figure indirectly compares the 870 /im visibilities and SEDs for these model structures with the 
observations. All of the model structures provide excellent fits to the broadband SED and Spitzer 
IRS spectrum. Individual SED model behaviors are indistinguishable, aside from small deviations 
near the far-infrared turnover where the measurement uncertainties are largest. However, the 
different structures and model types exhibit distinctive signatures in their resolved 870 /xm emission 
profiles. The similarity solution models with positive density gradients (7 > 0; Models sA and sB) 
substantially over-predict the emission on ~100-200 kA scales, while those with negative gradients 
(7 < 0; Models sD and sE) have visibility nulls at ~150kA that are clearly not commensurate 
with the data. With an intermediate 7 = 0, Model sC provides the best match to the continuum 
data for this model type. However, it and all other similarity solution models fail to reproduce the 
visibility oscillations beyond 200 kA. This "ringing" is a classic sign of a sharp edge in the radial 
emission profile, which is naturally produced with the power-law models described by Eq. (3). 
For those structures, steep density gradients {p = 1.0-1.5, Models pA and pB) produce too much 
emission on 100-200 kA baselines and shallow gradients {p = 0.0-0.5, Models pD and pE) generate 
visibility nulls that are not observed. However, an intermediate case with p = 0.75 (Model pC) is 
an excellent match to the data, with only a small (albeit significant) departure on 280-380 kA scales. 
That mismatch is likely related to the shape of the outer edge, although we have not pursued that 
speculation further. Of all the dust structures explored here. Model pC is clearly favored. 

A consistent feature of all the dust-based model structures is their compactness. The similarity 
solution models require characteristic radii of rc = 35-45 AU and the power-law models call for sharp 
outer edges at rg = 51-77 AU. That compact dust distribution was noted in §3, where we pointed 
out that all of the 870 /xm dust continuum emission is concentrated inside a radius of roughly 1" 
(~60AU; see Figure [Dj). Given the much larger radial extent of the CO J=3— 2 emission - out to 
radii of at least 4" (215 AU; see Figure [T]i) - the best-fit models face major problems reconciling 
the CO and dust observations. The moment maps in Figure [5] confirm this tension between the 
dust-based structure models and CO data, highlighting a clear CO-dust size discrepancy in nearly 
all cases. More direct comparisons of the CO channel maps with each model structure can be made 
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Fig. 5. — Moment maps of the CO J=3— 2 emission from the TW Hya disk and the various disk 
structure models compiled in Table [2j The leftmost panels show the SMA observations. The top 
panels make a direct comparison with the similarity solution models, and the bottom panels do the 
same for the power-law models with sharp edges. In all plots, contours mark the velocity-integrated 
CO intensities (0*'^ moment) at 0.4 Jy km s~^ (~3(t) intervals and the color scale corresponds to 
the intensity- weighted line velocities (l'^* moment). Only Model sA provides a good match to the 
observed CO emission; all others predict gas distributions that are too small relative to observations. 
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Fig. 6. — A direct comparison of the observed CO J=3— 2 channel maps with synthetic data 
from the similarity solution models in Table [2j The top block of panels show the model spectral 
visibilities synthesized in the same way as the data, and the bottom block displays the imaged 
residual visibilities. All maps have the same contour intervals as in Figure [Hi. The viewing 
geometry of the TW Hya disk is marked with a gray cross in each channel map. 
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Fig. 7. — Same as Figure [U but for the power-law models with sharp outer edges. 
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in Figures [6] and [71 which show both the models {top panels) and the imaged residual visibilities 
{bottom panels) synthesized in the same way as the data. Thanks to the freedom afforded by the 
parametric treatment of the gas temperatures, all models successfully reproduce the central CO 
emission core that is generated in the line wings. However, only Model sA has a sufficient amount 
of mass at large disk radii to account for the observations near the systemic velocity. Unfortunately, 
Model sA provides a poor match to the resolved 870 /im continuum emission profile. 

Before investigating potential physical explanations for the apparent CO-dust size discrepancy, 
it is important to evaluate the possibility that this is an artifact of some assumption in the modeling 
(aside from the underlying parametric models themselves; see §5). Admittedly, we adopted a relaxed 
approach for treating the relative vertical distributions of gas and dust in these models. Only two 
basic requirements were enforced for consistency: the gas and dust are in thermal equilibrium at the 
disk midplane, and the gas is distributed vertically according to hydrostatic pressure equilibrium. 
There are more sophisticated treatments of the gas in disk atmospheres that accurately account 
for the effects that X-ray heating, chemical diffe rentiation, and gas-dust temperature departures 
have on the vertical structure of the disk (e.g., 



Kamp fc Dullemondll2004l : iKamp et al.ll20ld : 



van Zadelhoff et al. 



Woitke et al 



200C 



2001 



Jonkheid et al. 2004 



Aresu et al 



2011 



A nearby 

disk like TW Hya would benefit fron i a more detaile d analysis with such models, particularly using 
multiple resolved molecular lines (see lThi et al.ll2010l . for a start). With a more physically motivated 
treatment of the vertical structure, we might infer a quantitatively different set of atmosphere 
properties that would affect the normalization and shape of the CO emission profile (currently 
controlled by {Ti, q} and other fixed parameters; see §4.2). However, the added complexity of 
those models would not change the fact that a compact density distribution is required to explain 
the dust emission. And since we assumed that the CO gas traces the dust in the radial direction, 
the size discrepancy would remain: these dust-based density profiles do not have enough material 
to be able to excite CO emission lines at large disk radii (see Figure [3]). 

To summarize, the fundamental conclusion from the radiative transfer modeling analysis of 
the SMA data is that the spatial distributions of the CO and millimeter-sized dust in the TW Hya 
disk are different: the dust is systematically inferred to be more compact than the CO gas. 



5. Discussion 

We used the SMA to observe the 870 fim continuum and CO J=3— 2 line emission from the 
disk around the nearby young star TW Hya. These data represent the most sensitive high spatial 
resolution (down to scales of 16 AU) probes of the millimeter-wave dust and molecular gas content 
in any circumstellar disk to date. Along with the SED, these observations were used in concert 
with continuum and line radiative transfer calculations in an effort to extract the disk density 
structure. We were unable to identify a consistent model structure that simultaneously accounts 
for the observed radial distributions of CO and dust. Assuming a radially constant grain size 
distribution and (vertically integrated) gas-to-dust mass ratio, the millimeter-sized dust structure 
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is significantly more compact than the CO. The resolved continuum emission profile demonstrates 
that the radial distribution of the millimeter-sized solids in the TW Hya disk has a relatively sharp 
outer edge near 60 AU, which is considerably smaller than the observed extent of the CO emission 
(out to at least 215 AU). 

The TW Hya disk structure has been studied extensively with millimeter-wave observations 
at lower angular resolution. In . a seri es of investigations focuse d almost exclusively on spectral line 



emission, 



Qi et al.l (|2004l . l2006l . l2008l ) and Hughes et al 



(120 llh relied on a slightly modified version 



of the physical structure model that was developed by ICalvet et al.l (j2002l ) to match the TW Hya 
SED. While that model has been used quite successfully to explain the molecular gas structure, there 
has always been som e tension between observations and its p redicted millimeter /radio continuum 



emission profile (see lOi et al.l l2004l : IWilner et al. 



2003 



20051 ). Given the modest quality of the 



previous continuum data and the fact that this structure model was not designed (or fitted) with 
access to resolved observations of any kind, the disagreement was understandably dismissed. As 
might be expected from its assumpti on of a "steady" a ccretion disk structure (i.e., with a large, 
positive surface density gradient), the lCalvet et al.l (j2002l ) model exhibits the same type of behavior 
as Models sA/B or pA/B: it agrees with the data on large scales (<100kA), but significantly over- 
predicts the amount of emission on smaller scales. The s a me is true for the parametric similarity 
solution models (where 7 « 1) explor ed bv iHughes et al.l (120081. l201ll). and w ould certainly apply 



to the analogous models developed bv lThi et all (j20ld ) and lGorti et al.l (j201ll ). All of this previous 



modeling work admirably reproduces the extended molecular line emission that is observed, and in 
many cases the SED and millimeter-wave continuum emission on large angular scales as well. It is 
only with the sensitive, high angular resolution data presented here that we recognize a problem: 
the radial distribution of the millimeter-sized dust grains is much more compact than for the CO. 

Unlike the thermal radiation at millimeter wavelengths, the optical and near- infrared light that 
scatters off small (<1 fim) grains in the TW Hya disk surface is detected out to large dis tances from 
the central star - at least ^ 4", comparable to what is inferred from CO spectral images (Krist et al 



2OOOI : iTrilling et al.ll200ll : IWeinberger et al.ll2002l : lApai et al.ll2004l : iRoberge et al.l l200,5l ) . So some 
dust traces the molecular gas, even if it is only a limited mass of small grains up in the disk 
atmosphere. However, these ex quisitely detailed scattered light images exhibit subtle structural 
complexities. iKrist et al.l (j2000l ) identified four distinct radial zones in the optical scattered light 
disk, with a prominent steepening of the brightness distribution just outside a radius_o£50AU_(~l''; 



their Zone 1/2 bounda r y). S imilar infrared behavior is noted in the studies by [Weinberger et al 



(2002) and 


Apai et al. 


( 


2004) 



80 AU (^1.0-1.5^^ ) range. Those results were confirmed in a comprehensive analysis of new data 



bv lRoberge et al.l (|2005l ). who also called attention to a color change at a similar radius as well as 
an azimuthal asymmetry out to a slightly larger distance from the star (~135AU). The physical 
origin of these scattered light features has been a mystery, although speculation centered around 
variations in the dust height (shadowing) and gradients in the dust scattering properties (either 
mineralogical or size-related) . But in light of our discovery of an abrupt drop in the millimeter- wave 
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continuum emission at the same location as these features, it is only natural to suspect that a more 
fundamental change occurs in the physical structure of the TW Hya dust disk near 60 AU. 

Perhaps the most straightforward explanation of the apparent CO-dust size discrepancy in- 
ferred from the SMA data is that we have used an incomplete description of the disk structure. 
As an example, consider a modification of either model type that incorporates an abrupt decrease 
in the surface densities (or millimeter-wave dust opacities) - not the dust-to-gas ratio - outside 
r PS 60 AU. If that drop in S (or Kmm) was not too large (perhaps a factor of ~100), there would still 
be enough disk material to produce bright emission from the optically thick CO lines and scattered 
light while also accounting for the sharp edge feature noted in the optically thin 870 fxm emission 
profile. That "substructure" in the outer dust disk might actually enhance the local gas-phase CO 
abundance, as ultraviolet radiation can penetrate deeper into the disk interior and photodesorb CO 
from the (small) reservoir of cold dust grains that remains at large radii (e.g.. iHersant et al.ll2009l ). 



Nevertheless, a physical origin for such a dramatic drop in the dust densities and/or the 
millimeter-wave dust opacities is not obvious. One possibility is that the disk has been perturbed by 
a long-period (as yet unseen) companion. If a faint object is embedded in the disk near the apparent 
edge of the 870 fim emission distribution, it might open a g ap that splits the disk into two distinct 
reservoirs and generate the warp asymmetry suggested by iRoberge et al.l (|2005l ). But, a narrow 
gap alone would not account for the SMA continuum observations. The millimeter- wave luminosity 
exterior to the gap would still need to be decreased, perhaps because the parti c les at those larger 
radii were preferentially unable to grow to millimeter sizes. IWeinberger et al.l (l2002l ) quote deep 
limits on i/-band point sources in the TW Hya d isk that sugge s t the re are no companions more 
massive than ~6 Mjup near 60 AU, acco rding to thelBaraffe et al.l (j2003l ) models (the corresponding 
mass limit is higher for the models of iMarlev et al.l 120071 ) . Certainly this kind of truncation or 
other forms of substructure could be invoked to explain the sharp radial edge in the SMA dust 
observations. But rather than engage in further speculation on the details, it should suffice to point 
out that the potential for substructure or other anomalies in the TW Hya disk can be tested with 
a substantial increase in resolution and sensitivity. Moreover, spectral imaging of optically thinner 
gas tracers (e.g., the CO isotopologues) would make for an ideal test of the origins of the apparent 
CO-dust size discrepancy. Fortunately, such observations will shortly be available as the Atacama 
Large Millimeter Array (ALMA) begins routine science operations. 

There is a compelling alternative explanation that has a more concrete physical motivation. 
In any protoplanetary disk, the thermal pressure of the gas is thought to cause it to orbit the 
star at slightly s ub-Keplerian rates, g enerating a small velocity difference relative to the particles 
embedded in it (IWeidenschillingj 119771 ) . Depending on their size, particles can experience a head- 
wind from this gas drag that decays their orbits and sends them spiraling in toward the central 
star. This radial drift of particles modifies the radial dust-to-gas mass ratio profile and introduces 
a pronounced spatial gradient in the particle si z e distribution - in es sence, it causes ( and Kmm to 
decrease with radius (IBrauer et al.l 120071 . 120081 : iBirnstiel et al.l l2009l). In the outer disk, t he drift 
rates are expected to be largest for millimeter-sized particles (e.g., iTakeuchi &: Linll2002l ). Since 
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thermal emission peaks at a wavelength comparable to the particle size, the drift process should 
then naturally produce a millimeter-wave er aission profile that is c onsiderably more compact than 
would be inferred from tracers of the gas ([Takeuchi &: LinI |2005| ) . Moreover, the dust particles 
that reflect light in the disk atmosphere are small enough to be dynamically coupled to the gas - 
therefore, scattered light images should be extended like the probes of the gas phase. 

From a qualitative perspective, our analysis of the CO and dust structures in the TW Hya disk 
are certainly consistent with a scenario where growth and radial drift have had an observable impact. 
However, it is unclear if realistic models of the growth and migration of solids in a disk like this 
can quantitative ly account for the deta ils. One particular challenge worth highlighting is related to 
timescales. The lTakeuchi &: LinI (|2005l ) models suggest that millimeter- wave dust emission should 
be strongly attenuated on a timescale of ~1 Myr without constant replenishment (presumably from 
growth and/or fragmentation). Given the advanced age of TW Hya (~8-20 Myr), this model would 
require either that replenishment shuts off after several Myr or that drift is inefficient at early 
times before becoming more important later in the disk evolution process. If this is indeed the 
mechanism responsible for our results, additional observations of disks at a range of ages could be 
used to help calibrate models of the long-term evolution of drift rates. One potential way to test 
this hypothesis relies on the particle size dependence of the radial drift rates. The theory implies 
that larger particles will end up with more centrally concentrated density distributions. At long and 
optically thin wavelengths, we expect that the size of the continuum emission region should be anti- 
correlated with wavelength - longer wavelengths imply more compact emission. In principle, this 
could be tested in the near future by combining high resolution ALMA and Expanded Very Large 
Array (EVLA) observations of the TW Hya dust disk that span the millimeter /radio spectrum. 

In many ways, TW Hya and its disk are unique and may not be representative of the bulk 
population of pre-main sequence stars and their circumstellar material. Nevertheless, it is worth 
considering the broader implications of our findings for this specific example - they may prove 
to be more generally applicable. It is possible that the CO-dust size discrepancy found here is 
present in most other millimeter-wave disk observations, but it would likely be difficult to identify 
with current sensitivity and resolution limitations. If that feature is common and its underlying 
cause is a drop in the dust-to-gas ratio in the outer disk, there would be serious consequences for 
disk mass estimates based on dust continuum measurements. There could be large and hidden 
mass reservoirs of molecular gas in the outer reaches of protoplanetary disks, implying that disk 
masses might be substantially under-estimated. If true, gas densities at large disk radii may be 
higher than typically assumed, with profound implications for fa c ilitating gi ant planet formation 
by gravitational instability (e.g.. iBolevi l2009l : iKratter et al.l 120101 : IBossI I2OIII ). However, if radial 
drift is the key process responsible for the apparent CO-dust size discrepancy, it is also possible 
that any "primordial" dust-to-gas ratio integrated over the entire disk is preserved. This would 
be the case if the millimeter-sized dust originally present at large radii had its inward migration 
halted before it was accreted onto the central star. In that scenario, the total disk mass estimates 
from millimeter-wave luminosities would still be relatively accurate (assuming a proper model for 
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the dust opacities that considers the simultaneous particle size evolution is available), although 
the densities in the outer disk would remain uncertain without measurements of optically thin gas 
emission lines. 

Independent of the apparent CO-dust size discrepancy, our finding that the millimeter-wave 
continuum emission from the TW Hya disk is so sharply truncated comes as a surprise. Modern 
interferometric datasets generally do not have sufficient sensi tivity to differentiate between dust 



density models with sharp edges or smooth tapers (e.g., see llsella et al.l l20ld : iGuilloteau et al 



201l|). Given that ambiguity, it is possible that the edge feature we have identified in the 870 /im 
emission profile of the TW Hya disk could be relatively common. Moreover, it is tempting to 
associate this r ~ 60 AU edge in the distribution of larger solid particles in the TW Hya disk with 
the similarly abrupt truncation of the cla, ssical Kuiper Belt at r ~ 40-50 AU in our solar system 



(jTrujillo et al 



2001 



Gladman et al 



2001 



If this behavior ends up being a generic feature in 
protoplanetary disks, it may signify an important diagnostic of the radial migration of disk solids 
an d provide new insight s into the structural origins and evolution of the outer solar system (e.g.. 



see 



Kenvon fc LmJll999l : iLevison fc Morbidellill200.'j ). 



Further speculation on the generality of the features identified in the TW Hya disk is unnec- 
essary. With the recent start of ALMA science operations, the quality of the data presented here 
will be matched (and exceeded) routinely for large samples, and the basic trends of disk properties 
like those probed here will be clarified. If the TW Hya disk is not anomalous, it is clear that the 
general methods used to interpret observations of dust in disks will need to be modified to focus 
less on their bulk density structures and more on the dynamical evolution of their solid contents. 



6. Summary 

We have presented sensitive, high resolution (0'.'3 = 16 AU) SMA observations of the 870 //m 
continuum and CO J=3— 2 line emission from the disk around the nearby young star TW Hya. 
Based on two different parametric formulations for the disk densities, we used radiative transfer 
calculations to compare the predicted radial structures of the dust and CO gas in the TW Hya disk 
with these SMA observations and ancillary measurements of the spectral energy distribution. The 
key conclusions from this modeling analysis of these high-quality data include: 

1. Under the assumption that the dust-to-gas surface density ratio is constant with radius, we 
were not able to find any model structure that can simultaneously reproduce the resolved 
brightness profiles of the 870 fxm continuum and CO line emission. We have identified a clear 
CO-dust size discrepancy that is present regardless of whether the assumed surface density 
profile has a sharp outer edge or a smooth (exponential) taper at large radii. 



2. The radial distribution of millimeter-sized dust grains in the TW Hya disk is substantially 
more compact than its CO gas reservoir. The 870 /xm dust emission has a sharp outer edge 
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near 60 AU, while the CO emission (and optical/infrared scattered hght from smah grains 
that are dynamicahy coupled to the gas) extends to a radius of at least 215 AU. 

3. The observationally inferred CO-dust size discrepancy could potentially be explained with a 
more complex dust density profile that exhibits a sudden decrease by a large factor near a 
radius of 60 AU. That "break" in the density profile might be consistent with a tidal pertur- 
bation by a long-period (unseen) companion, although the dust at larger radii would then 
also have to be preferentially depleted of millimeter-sized grains. 

4. Alternatively, the observations might have uncovered some preliminary evidence for a key 
evolutionary mechanism related to the planet formation process: the growth and inward 
migration of disk solids. The radial drift of millimeter-sized particles is expected to naturally 
concentrate long- wavelength thermal emission near the star relative to tracers of the molecular 
gas reservoir. However, a more detailed exploration of disk evolution models is needed to verify 
if the observations of the TW Hya disk are quantitatively consistent with this scenario. 



Wc thank Mark Gurwcll for his assistance with the multi-epoch calibration of the SMA data, 
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